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Currently laser cooling schemes are fundamentally based on the weak couphng regime. This 
requirement sets the trap frequency as an upper bound to the cooling rate. In this work we present 
a numerical study that shows the feasibility of cooling in the strong coupling regime which then 
allows cooling rates that are faster than the trap frequency with state of the art experimental 
parameters. The scheme we present can work for trapped atoms or ions as well as mechanical 
oscillators. It can also cool medium size ions chains close to the ground state. 



Introduction — Laser cooling is the main tool that en- 
ables the exploration of low temperature phenomena in 
atomic physics, at nano and micro systems as well as the 
rapidly developing field of quantum technologies. Fol- 
lowing the original ideas of Doppler cooling 1, 2], laser 
cooling has taken a central role in the physics of cold 
atoms and the number of ideas and their sophistication 
is growing continuously. In the quest to propose cooling 
schemes that can cool to lower and lower temperatures 
at ever increasing rates the complexity and the efficiency 
of laser cooling has progressed a long way. 

Two level systems(TLS) admit the Doppler cooling 
limit. This limit can be overcome by the Sisyphus mech- 
anism either for free0-Q or trapped particles^]. The 
recoil limit can be broken by higher level systems by Ve- 
locity Selective Coherent Population(VSCP)f7'| trapping 
or by Raman cooling [8] for free particles and in their ana- 
log for bound particles, Dark State cooling 9j and Raman 
side band cooling Cooling schemes for trapped parti- 
cles can be applied for mechanical oscillators in the setup 
described in [ll|, [12 1 . 

However, all such cooling schemes rely on a weak cou- 
pling between the internal degrees of freedom (dof) and 
the ion's external motional state. This results in a clear 
separation of time scales between the two, allowing adia- 
batic elimination of the faster of the two dof-s. Unfortu- 
nately due to the weak coupling nature, the trapping fre- 
quency necessarily becomes the upper limit to the cooling 
rate and thus the achieved cooling rates are a few orders 
of magnitude lower than the trap frequency. 

We here present a fundamentally novel way of cool- 
ing a system of trapped particles, achieving rates that 
are faster than the trapping frequency, and which grows 
with the laser intensity with no fundamental upper limit. 
We demonstrate by means of numerics the proposed prin- 
ciple on a linear ion trap and show in detail how, using 
optimized sequences of coupling pulses, rapid cooling can 
be achieved. The proposed method is robust to fluctu- 
ations of the laser and may be adapted to a very wide 
range of systems, this can be achieved, for neutral atoms 
Q , trapped ions at low temperature j^] , and mechanical 



oscillators [ll|, |l2| for which the method can break the 
temperature limit imposed by the finite Q factor. 

The superfast cooling concept — Consider a trapped 
system which is coupled to a TLS. The system is in an 
almost-harmonic potential, with engineered coupling be- 
tween the external dof and the internal ones. For a stand- 
ing wave configuration when the ion or atom is at the 
node we obtain in leading order the following Hamilto- 
nian: 
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where 77 is the Lamb-Dicke parameter, v is the trapping 
frequency, is the effective Rabi coupling between the 
two internal levels, a is the phononic annihilation oper- 
ators and <7g is the pseudo spin operator in the direc- 
tion 9 in the x, y plane. To cool the system, we would 
like to transfer energy from the motional to the internal 
dof., which can be dissipativly reinitialized via optical 
pumping. In other words, we wish to implement the red 
sideband-cooling operator, aa^ -|-aV~, \2^^2^]. In side- 
band cooling this term is generated from the Xgx term 
after a rotating wave approximation (RWA), which re- 
quires a time which is longer than the trapping frequency, 
where X is the position operator. 

A natural question to ask is whether we can cool, i.e. 
create the sideband-like term, in times which are faster 
than the trap frequency. A positive approach is by nu- 
merical optimization, but this requires an important in- 
put which is the initial 'point', i.e., a guess of the cooling 
cycle. One may try to initiate the optimization with 
the schemes valid for the weak coupling regime such as 
sideband cooling or the dark state cooling schemes and 
extend them to strong couplings. This approach proves 
highly inefficient. In our work we approach the prob- 
lem from the opposite direction: we choose as our initial 
guess a cooling scheme in the impulsive limit (infinitely 
fast), where we could reach the ground state infinitely 
fast if we neglect the higher order terms in the Hamilto- 
nian. Starting from this point we optimize numerically 
to adapt finite couplings attainable in the lab. 

We start by discussing the idealized case to gain in- 
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sight into the optimization procedure. Coohng at the 
strong-couphng hmit is simple and extremely fast. In 
this limit there are various 'points' from which we can 
start the numerical optimization. We choose to start 
with the following 'point'. Consider the following ar- 
gument which elucidates the underlying intuition of our 
work: The side band cooling term can be written as 

aa^ + a^ffj- = (^Xa^ - :;^P(Ty^ ■ We already have 

a Xax coupling available to us ([ij . Thus if we can create 
a Pay term, we can generate the red sideband Hamil- 
tonian using the Trotter decomposition [17[, neglecting 
some constants: (^e'^Xa^dt^tnPaydty ^ ^i{xa^-p<yy)Q 

when dt ^ and flndt = Q; i.e., multiple short pulses of 
Xax and Pay will create the needed Hamiltonian. There- 
fore, if we had infinitely strong lasers and the Pay in- 
teraction, we would be able to cool instantaneously, at 

e = TT. 

Pay term — Following the insight first described in 
(l8| we derive Pay as an effective Hamiltonian using Xay 
pulses which give the ion momentum, a period of free 
evolution in which the ion translates and finally a —Xay 
pulse, imparting equal and opposite momentum, stop- 
ping the translation. Furthermore, in the strong coupling 
regime the coupling is much larger than the phonon en- 
ergy, ryfi ^ v, allowing us to ignore the free evolution 
for the duration of the pulses. Setting (5 = yields a 
solvable Hamiltonian which enables us to get an exact 
result. Formally, splitting ([1} into Hjree = [a^a + ^) 
and Hpuise = hilrj {a) + a) ay, setting 5 to zero and mak- 
ing use of the Baker-Campbell-Hausdorff formula (l9l - [2lj | 
we get a closed-form expression 
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which includes the desired Pay operator. We can now 
combine a Xax pulse with a Xay-wait-Xay sequence im- 
plementing a Pay pulse and generate a red-sideband op- 
erator. Note the we have a large phase term firf'vVL^t ft^ 
that we will have to deal with. Initial numeric opti- 
mization assumed that idealized P pulses are possible, 
with relatively small density-matrix cut-off sizes, which 
together allowed a much-accelerated computation proce- 
dure, and then further work adapted these sequences to 
the realities of achievable pulse strengths and validated 
the cooling processes with much larger d.m.-s. As ex- 
pected, the idealized pulses resulted in even better per- 
formance than is presented in this letter. 

Numerical optimization of realistic scenarios — Mov- 
ing closer to a true model of an ion-trap, we need to move 
away from the impulsive limit and re-introduce free evo- 
lution during pulsing, as well as account for residual de- 
tuning and higher-order elements in the Lamb-Dicke ap- 
proximation. In this case the Hamiltonian is not solvable 
anymore and the BCH form employed for the Pay denii- 
pulse does not apply. Thus a correction should be added 



to Eq. [21 The free evolution during the pulse changes 
the commutation relation between the two pulses as the 
exponent in the second pulse is not the negative form of 
the exponent in the first pulse - the free evolution does 
not change sign. This gives a crucial correction to Eql^l 
which has to be taken into account. Here we have a 
major problem, as this correction cannot be taken into 
account in a purely numerical way since this will force 
us to introduce numerically the operation of the strong 
pulse which will require us to use a truncation at large 
phonon numbers. 

Moreover, the Pay demi-pulse generation is deeply 



connected to the 



X,P 



commutation relation, which 



breaks once one tries to use finite-sized matrices in nu- 
merical calculations, and the size cut-off introduces sig- 
nificant additional inaccuracies when trying to make use 
of these commutation relations via 3 matrix exponents. 
We must therefore compute the equivalent to the pulse- 
wait-pulse sequence analytically, merging the 3 exponents 
into one, and only then can we safely chain operations 
into a complete cooling sequence, which is less sensitive 
to numerical inaccuracies. We therefore employ the ex- 
panded BCH form as appearing in f22|. 

Unlike the BCH form used to derive the effective P- 
pulse in the infinite limit, the BCH series does not ter- 
minate after a finite number of elements, and has to be 
computed to high order to achieve accurate results. To 
achieve a workable, reasonably-compact form for such a 
series requires the use of computerized algebra which is 
capable of making use of commutation-relations to sim- 
plify expressions. This has been implemented in Mathe- 
matica. This avoids errors due to [x,p\ numerical failure. 

Results — As the Pay interaction shifts the loca- 
tion of the ion it cannot be performed instantaneously. 
Moreover, the true impulsive limit H free'^Hpuise is not 
currently accessible in the lab and therefore we can- 
not completely ignore free evolution while pulsing. As 
a result, while the above approach provides us with a 
framework and a starting point, we must employ quan- 
tum optimal control techniques, the details of which are 
achieved through numerical optimization, to apply the 
above methods for finite pulse lengths and finite coupling 
strength. We combined alternative steps of simulated an- 
nealing, which excels at extricating the search out of local 
minima but is slow at convergence to the optimum, with 
BFGS gradient following search, which excels at fast con- 
vergence but has an observed tendency of getting stuck 
at local minima. 

Numerical studies for this work have been per- 
formed using QLibTs*!, a Matlab package for quantum- 
information and quantum-optics calculations. The pulses 
were optimized to give the lowest possible average num- 
ber of phonons after cooling. 

We optimized cooling cycles for a given initial temper- 
ature. The optimization was performed in the impulsive 
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Initial energy [hi/] 
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l<inal energy [hu\ 


0.4 


1.27 


1.95 


Final energy after 25 cycles 


0.02 


0.10 


0.22 


Cycle duration 


4 4 


2.7 — 

71 


0.8 ^ 


No. pusles per cycle 


180 


90 


90 


No. of sequences 


10 


10 


10 



TABLE I. Cooling cycle performance. We define a cooling 
sequence as a series of alternating Xa^j-pulses and Pay-demi- 
pulses of varying lengths, followed by a reinitialization of the 
ion's internal dof. A cooling cycle is comprised of several 
cooling sequences, which are generally non-identical. 
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FIG. 2. Repeated applications of the sample cycles. 
Inset: X and P pulse lengths for sample cooling cycle C. Neg- 
ative pulse lengths indicate the spin is to be oriented in the 
negative x or y directions. Each block of 6 bars (3 X and 3 P 
pulses) represent a single cooling sequence. These sequences 
are interspaced by reinitialization of the ion's internal state 
to make up a complete cooling cycle. 



FIG. 1. Initial and final energy (above ground, in units of 
hi/) for a single application of the cooling cycle and for 25 
applications. 



limit (rjD, ^ v), and applied for realistic parameter val- 
ues, where we may not neglect free evolution while puls- 
ing. The sequences presented are just examples of what 
is achievable, and by no means are they to be considered 
canonical in any way. 

All computations below are done for the following 
physical settings, which are achievable in the lab: *'^Ca^ , 
V = IMHz • 27r, = IQQMHz ■ 27r, Xiaser = 730nm, giv- 
ing a Lamb-Dicke value oirj = 0.31, IJ]. Table 1 details 3 
sample cycles and their cooling performance. Note that 
the coupling terms in the propagator are proportional 



to tpT]^ ^^^{IXag and rjiftpflvPay and therefore total 

cooling time scales as I/a/TJ since tp x tf scales like l/fl. 

We made the optimization in the impulsive limit, and 
then applied the deduced sequence using the full Hamil- 
tonian. The results are shown in Fig. [1] Starting from a 
thermal state of varying temperatures, Fig[T^ shows the 
final energy for a single application of a cooling cycle; 
Fig lb shows the result of 25 applications of the same 
cycle. The latter cannot be deduced directly from Fig la 
because the state is no longer thermal after one or more 
cooling cycles. Note how the cycles give a good cool- 
ing for a wide range of initial temperatures, even though 
optimization was carried out for a specific initial temper- 
ature 

In Fig. [5] we show how repeated applications of the 
three sample cycles continues to lower the energy, until 
a steady state, specific to each sequence is achieved. An 



example timing of sequences within a cycle is shown in 
fig. [2] inset. 

Robustness — The robustness of the cooling pulses is 
extremely important for experimental realization. In or- 
der to analyze the robustness we have simulated the op- 
eration of the pulses under noisy conditions of the lasers. 
The results indicate that the cooling pulses are extremely 
robust even though they were not optimized to be so. In 
Fig. [3] we study the robustness of cycle C, repeated 25 
times, to Gaussian noise in the pulse timings. We assume 
that instead of the prescribed pulse time t, we implement 
t — 5- (1 + e) t, with e being drawn from a Gaussian dis- 
tribution with a varying standard deviation. A similar 
noise pattern was applied to the laser power, il. The 
superfast cooling exhibits two slightly different sensitiv- 
ities to noise: for the short time-frame which 
interferes with the commutation relations allowing the 
creation of the effective P pulse, sensitivity is somewhat 
higher, while for noise occurring on longer timescales, 
i.e. the duration of the various pulses comprising the 
cooling sequence, Superfast cooling is extremely robust. 
The plots show the mean final phonon count over 500 cy- 
cles. A graceful degradation in cooling performance can 
be observed. 

Note that the state of the system after cooling is not 
a thermal state. This is due partly to purely numeri- 
cal issues stemming from the need to exponentiate what 
are essentially infinite-sized Hamiltonians, and partly a 
physical feature of the specific cooling sequences pre- 
sented. The latter results from the necessity to make use 
of relatively small density matrices in initial optimization 
stages, which resulted in effective cooling operators that 
do not cool high harmonic modes (above 30) as well as 



4 




123456789 10 



Relative Gaussian noise [%] 

FIG. 3. Robustness to noise of cooling cycle C, applied 25 
times. The x - axis represents the relative noise and the y- 
axis the final phonon population. As can be seen, superfast 
cooling is proven to be quite robust, both at long (u) and 
short (<C I/) timescales. 
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Number of ions 

FIG. 4. The average population of the center of mass mode 
as a function of the number of ions, o is for a regular chain 
and o is for pinned ions. 



they cool lower modes. While this difference is negligible 
before cooling, when going significantly below 1 phonon, 
the remaining energy in the high modes becomes a not- 
insignificant part of the overall energy. We believe that 
with higher computational power sequences resulting is 
more thermal-like final states are achievable. 

Ion Chain — - By applying this cooling method to 
a large chain a product state of ground states can be 
reached, i.e., we will cool to the |0) |0) .... |0) state in the 
local basis. This state is not far from the global ground 
state. In fig. 2] we show the number of phonons in the 
center of mass mode for a regular ion trap and for a trap 
in which all the ions are equidistantly pinned. In the 
latter case a saturation value of around 0.2 phonons is 
reached for the center of mass mode while for the regular 
trap there is no staturation but the population is still low 
for moderate chains. 

Standing wave — - There are a few ways to gener- 
ate the standing wave Hamiltonian. We are not re- 
ally interested in a standing wave but in a pulse of 
a standing wave. We can emulate this pulse by run- 
ning wave pulses by using the two following pulses: 
ifl {ux + rjx) , —iO. {ux — rjx). These pulses can be cre- 
ated by passing a pulse through a beam splitter and 
bringing it through two different directions, thus con- 
trolling the phase and sign of the Lamb Dicke param- 
eters independently. The unitary transformation which 
describes this is: e*"'"-e*''"''-="e-*"'"-e*''"''-^ = ^t2nna^x ^ 
which creates the unitary operation that we need. More- 
over, since we are interested in fast cooling rates we also 
are ensured that the standing wave will not drift. The 
creation of the standing wave can be done as in[26']. In 
the case of mechanical oscillators this trick would not 



work since there is no temporal control of the effective 
Lamb Dicke parameter, but if the coupling is with Ra- 
man transitions a different intermediate level can be used 
and thus the effective Lamb Dicke will be different while 
the Rabi frequency may be the same and thus cancel the 
first pulse. 

So far we have assumed the ion can be modeled by a 
non-dissipative two-level system with effective Rabi fre- 
quency ri. This is a simplification. Some experimental 
setups use a 3-level system with the effective Rabi cycle 
achieved by way of a highly detuned Raman transition 
off a highly-dissipative level [23,-.25j. As a result, the true 
laser coupling of \QQGHz is reduced to an effective co- 
herent lOMHz transition. However, unlike quantum in- 
formation processing, superfast cooling does not require 
a fully coherent transition. We could reduce the detuning 
and achieve a higher Rabi frequency. At = IGHz, for 
example, the cooling process speeds up, requiring only 

0. 4, 0.3 and 0.2^21 for cycles A, B and C, respectively. 

Nanomechanics — This scheme may have increasing 
importance in the nonmechanical and hydromechanical 
world since there the finite Q value limits the final tem- 
perature. One would use the setup described in 11. 1^: 

1. e. couple a TLS to the oscillator. In this case the ef- 
fective Lamb Dicke parameter is on the order of 7y = 
0.03([l2|) and can be much higher for the scheme which 
is generated with magnetic gradients ([lH). In both of 
these schemes the strong coupling regime can be reached 
by strong lasers or even microwaves. The method cur- 
rently in use 27|, |28| , assisted by a cavity, cannot beat 
the limit of the trap frequency due to the Gaussian na- 
ture of the scheme. 

Conclusion — To conclude, we have introduced the 
optimized-control approach to cooling trapped ions, by 
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optimizing over the realizable coupling between the 
phonons and the ion's internal levels, Xag, and periods 
of free evolution. Surprisingly, even this meager palette 
is sufficient to implement a red sideband-like operation 
and drive energy from the phonons to the ions faster than 
the trapping frequency. Furthermore, the principles de- 
scribed above are applicable to a very wide range of sys- 
tem and robust enough to be implementable in ongoing 
experiments. 
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